Chaos and its quantization in dynamical Jahn- Teller systems 
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We investigate the E g igi e g Jahn- Teller system for the purpose to reveal the nature of quantum 
chaos in crystals. This system simulates the interaction between the nuclear vibrational modes and 
the electronic motion in non-Kramers doublets for multiplets of transition-metal ions. Inclusion of 
the anharmonic potential due to the trigonal symmetry in crystals makes the system nonintegrable 
and chaotic. Besides the quantal analysis of the transition from Poisson to Wigner level statistics 
with increasing the strength of anharmonicity, we study the effect of chaos on the electronic orbital 
angular momentum and explore the magnetic g-factor as a function of the system's energy. The 
regular oscillation of this factor changes to a rapidly-decaying irregular oscillation by increasing the 
anharmonicity (chaoticity). 

PACS numbers: 05.45.Mt,71.70.Ej,82.90.+j,03.65.-w,31.30.Gs 
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I. INTRODUCTION 

Recently the study on quantization of classically 
chaotic Hamiltonian systems has received a wide atten- 
tion. An accumulation of numerical and experimental 
data indicates Wigner-type level statistics, wavefunction 
scars and other characteristic featurespIQ. 

In addition to toy models like a kicked rotator, Henon- 
Heiles system, etc., some realistic systems like a hy- 
drogen atom in a magnetic field and micro-wave cavi- 
ties are also being investigated^, 0, 0- Quantum me- 
chanics of chaotic systems also suggests insight beyond 
a simple quantal manifestation of chaos 0, [3. There- 
fore it is crucial to have more and more experimentally- 
accessible quantum systems which exhibit chaos in its 
classical treatment. 

In this paper we choose the Jahn- Teller system simu- 
lating transition-metal ions embedded in the host crys- 
tals such as III-V semiconductors and halides crystals. 
Among them we consider E g ® e g model associated 
with the irreducible representation for the cubic sym- 
metry group, namely, the two-dimensional (2-d) lattice- 
vibration modes e g linearly coupled to doubly-degenerate 
electronic states E g 6J. This system has an adiabatic 
doubly-fold lattice potential with the conical intersection 
of the potential surfaces .whose geometric phase was one 
of topics some time ago0,|^. The lattice potential here 
can be harmonic or anharmonic. From the classical dy- 
namical viewpoint in the adiabatic limit, as shown be- 
low, the system with the 2-d harmonic potential is inte- 
grable, leading to regular motions, and on adding the an- 



harmonic term, it becomes nonintegrable and chaotic@. 
A systematic investigation of the quantal counterpart of 
classical chaos in these systems is desirable. Furthermore, 
since the model is a representative for paramagnetic ions, 
it is experimentally important to see the effect of chaos 
on the magnetic g- factor. This factor is an expectation 
value for electronic orbital angular momentum and mea- 
sures a degree of level splitting of highly excited states 
induced by the lattice-electron interaction. The oscillat- 
ing structure in energy dependence of (/-factor is expected 
to reflect the feature of the underlying classical dynamics. 

The organization of the paper is as follows: In Section 
II a model for Jahn- Teller E g e g system is proposed. 
Section III deals with the classical analyses of the model. 
Both the systems with and without anharmonic terms are 
examined. Section IV presents a quantization of the sys- 
tem together with level statistics. Section V is concerned 
with a proposal of the experiment to verify the quantum 
signature of chaos in the dynamical Jahn- Teller system. 
Some novel feature of ^-factor is explored there. The final 
Section is devoted to summary and discussions. 
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II. DYNAMICAL JAHN-TELLER SYSTEM 



We investigate the electronic states of degenerate E g 
orbitals of d-levels in transition-metal ions coupled with 
2-d vibrational modes e g expressed by coordinates Qi 
and Q'2- The E g ® e g model is the typical system show- 
ing dynamic Jahn- Teller effects (DJTE), which has been 
discussed in the field of magnetism for transition-metal 
ions 0, 0| . The Hamiltonian matrix H for this system is 



2 



expressed as 



H 



h 2 d 2 

2 [ 8Q 2 



dQ 2 2 ]I 



Qi 



(i) 



where I is the 2x2 unit matrix and V(Q) is a potential 
energy. The nuclear mass is set to be unity. The second 
term of is the so-called Jahn- Teller interaction, Hj-t 
with k the coupling parameter between electronic states 
and vibrational modes. Bases for electronic orbitals E g 
lying behind are (r\u) — u(r) = 3z 2 



v(r) 
by 



2 and (r|«) = 

-y 2 . When V(Q) is a harmonic potential given 



without anharmonicity. We explore the effect of chaos 
on the g-factor in the system with the anharmonicity. 
The dynamical Jahn- Tell er sy stem was also studied by 
Bulgac and Kusnezov^], Il8l IT^. |2C| in a system with 
the three-dimensional harmonic potential. However, we 
should note that the dimensionality of lattice-vibration 
modes characterized by the irreducible representation is 
two and not three according to the theory of a point- 
symmetry group applied to real crystals and that our 
model is a better reflection of the real crystaljfj. 



III. QUASI-CLASSICAL DYNAMICS AND 
CHAOS 



V (Q) = ~u 2 (Ql + Q 2 2 ), 



(2) 



the corresponding adiabatic potential for has the 
Mexican-hut shape in Fig^ where Qi = pcos9 and 
Q 2 = psin9. The potential minima lie at p = p = k/uj 2 
with an arbitrary value of 9. Namely the minima are in- 
finitely degenerate. The energy for the minima is k 2 /2uj 2 . 
Vibronic levels for the quantum Hamiltonian were 
discussed in numerical calculations using small dimen- 
sional Hamiltonian matrices [ijj. Recently this model 
has been investigated from a viewpoint of the geomet- 
ric On the other hand, the effect of the 
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In the first place, we shall analyze the quasi-classical 
counterpart of Hamiltonian Q , which is given by 



H=-(P 2 + P 2 ) + V(Q) + k(Q-a), 



(4) 



where the first term is a kinetic energy for classical vi- 
brational modes with coordinates Q = (Qi, Q 2 ) and the 
second one is the harmonic and/or anharmonic poten- 
tial; the third one is the quasi-classical form for the 
Jahn- Teller interaction where tr is Pauli matrices er = 
(cr x ,a y ,a z ). Noting that a space is independent of the 
real space, we choose a x ,a y and a z corresponding to 
02,0-3 and <7i, respectively. It is convenient to represent 
the quantum state with use of the density matrix p 



P = 

where r = Tr (pa) = 
Using the potential 



1/ 1 + z x — iy 
2 I x + iy 1 — z 



(5) 



(x,y,z) = (tx,z) is a real vector. 



V r (Q) = U (Q) + U 4 (Q) 



(6) 



with Vo and Va in J2| and the equations of motion 
derived from are 



FI G. 1: A diabatic potential of Mexican-hut shape, p = 
\J Q\ The potential has the degenerate minimum. 

trigonal fields expressed as the anharmonic term 

V A (Q) = -(b/3)(Q 3 1 -3Q 1 Q 2 2 ) (3) 

was also analysed as to some low- lying levels 0, • I n 
short, O'Brien investigated the system (JIJ with poten- 
tial V(Q) = Vb(Q) + Va(Q) in the low-energy approx- 
imation that p is fixed to pg. However, we numerically 
calculate eigenvalues and eigenvectors without having re- 
course to such an approximation. We derive level spac- 
ing distributions to see the effect of chaos on quantum 
systems 0, 0, EJ 0, El Furthermore we investi- 

gate the quantal (7-factor, whose oscillating structure was 
shown three decades ago by Washimiya in the system 



dQ 
~~dt 
dP 

dt 



dr , _ 

— = fcQ x r. 

dt 



dV(Q) 
dQ 



kr 



(7a) 
(7b) 
(7c) 



Equation fTcjl is nothing but the Schrodinger equation 
ip = [H, p] , from which we find the constant of motion 
I r |= 1. In the study of quasi-classical dynamics, our in- 
terest lies in qualitative comparison between the systems 
with and without the anharmonic term, and therefore we 
confine to the adiabatic limit 4| = 0, that is, r = (rj_, zq) 
with r± 2 = 1 — z 2 . Further, we find from J7cJ), 



dz 
'dt* 



fcQ xri = 



3 




which is satisfied only when rj_ || Q. Thus the adiabatic 
limit is equivalent to 



r, 

r_L = VI ~ '- 



Q' 



where Q = |Q|. Thanks to ©, (EEJ) reduces to 

dP _ dV(Q) ~Q 
~dt ~ dQ Q 



(8) 



(9) 



with the renormalized coupling k — fcVl — zq 1 . Conse- 
quently, the classical equation of motion for the present 
model can be expressed only by Q, P and consists of a set 
of (|7a|l and ©. This set has the first integral of motion 
or the total energy 



P 2 ~Q 2 



(10) 



The following analysis depends on the type of the po- 
tential V^Q). Firstly, we investigate the system with 
the harmonic potential only, i.e., ^(Q) = Vo(Q). In this 
case, in addition to the total energy (|1U[) . we have another 
constant of motion, i.e., the orbital angular momentum 



Jz 



(Q 



(ii) 



The number of constants of motion agrees with the de- 
grees of freedom(two) . Therefore the system is integrable 
HI HH HHl - showing only regular motions. 



Then, we investigate the system with the anharmonic 
potential V(Q) = Vb(Q) + V A (Q). The trigonal field 
on the 2-d plane (Qu Q2) is invariant only to operations 
of the cubic group[6|- Owing to this breaking of con- 
tinuous circular symmetry, the angular momentum J z in 
(|llfl is not a constant of motion, which makes the system 
nonintegrable. It should be noted: the lattice system 
without coupling with the electronic degree of freedom 
is identical to the Henon-Heiles system whose dynami- 
cal features have been intensively studied in a context of 
chaos 

theory[! E3. 



The present system has two control parameters, i.e., 
the coupling constant k between electronic and vibra- 
tional degrees of freedom and the nonlinearity parameter 
b responsible for the trigonal field. However, the simple 
scaling below lets them merge to a relevant single pa- 
rameter. Let the coordinates {Qi,Q2) be transformed 



to (qi,q 2 ) through J ~Q\ 



\Qi = Q2- The total 
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energy is then written as 

(12)' 

Next, define the scaled time r = (fc6) 1 ' 4 t and the scaled 
momentum p = Finally, by scaling the energy and 
the angular frequency as 

e = Eb 1 '^' 2 (13a) 

n 2 = J 2 l- 1 ' 2 b- 1 ' 2 , (i3b) 

the total energy is expressed as 

e=\v 2 + + \ V^ 2 + | - 9l <7 2 2 . (14) 

Eliminating k between (|13a|l and (|13b|l . we find e — 
Eb 2 , which tells that the enhancement of the an- 
harmonic term is equivalent to the increase of the scaled 
energy e under a fixed value of — . In our numerics be- 
low, we choose § = 1. Figure [21 shows Poincare sur- 
face of sections for energies e = 0.01,0.4 and 0.7. While 
in Fig|2fa) we find almost all trajectories to be regular, 
FigGIb) shows that Kolmogorov-Arnold-Moser (KAM) 
trajectories begin to collapse and some trajectories be- 
come chaotic. Finally in Fig|2fc) chaotic trajectories 
dominate almost all phase space. These results imply 
that the domain of chaos expands gradually with increas- 
ing the scaled energy or the anharmonicity under a fixed 
value of the scaled angular frequency. 

We have also solved Eq.(7) without use of the adia- 
batic approximation. Then, the degree of freedom be- 
comes three. In case that k and b are less than unity, the 
relative fraction of chaos in phase space is quite small: 
We call this behavior "partial chaos" . However, if both 
k and b are as large as 5, the global chaos as seen in 
Fig^c) appears. In case of b = 0, we find no indica- 
tion of irregularity. We shall now proceed to investigate 
a quantal counterpart of these classical features. 
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FIG. 3: (a) Histograms of level spacing distribution 
for k — 4, b = 0(without trigonal field). Curves of 
Poisson distribution and GOE are also drawn by 
open rectangles and crosses, respectively, (b) Inte- 
grated level spacing distribution with use of data 
in Fig.3(a). 



where n = 1,2,..., and m = n — 1, n — 3, . . . , — n + 
1- F n \m\{p) is the confluent hypergeometric function |24|. 
Thus, the basis wavefunctions for the present model are 
given by 

®n,r n = U ±( r )<l>n,m{p,0)- (17) 



IV. QUANTUM SYSTEMS 

In this Section we investigate the quantal manifesta- 
tion of chaos and of regular motions without resorting 
to the adiabatic limit. The construction of the energy 
matrix is as follows(see ^l] and [24|): The basis wave- 
function is described as the products of electronic wave- 
functions u±(r) and vibrational ones </>(Q). The former 
is given in terms of u(r) and v(r) below (JIJ as 

u±(r) = -^(u(r)±*i;(r)), (15) 

while the latter is the eigenstate with the eigenvalue 
Enm = nJUd for the 2-d harmonic oscillator: 

<Pn,m(P,8)=Fn\ m \(p)e imd , (16) 



Corresponding to QJ, the interaction matrix Hj~t is 
expressed as 

Hj_ T = V u (r)Q 1 + V v (r)Q 2 

= ^=[V u -(r)e ie -V u+ (r)e- ie ] (18) 

with the matrix elements of V u ± given by 

(«±|K+(r)|«±) = («±|K-(r)|«±) = 
K|K,±(r)|«±) =TV2k. (19) 

As for the lattice- vibration factor pe ±l9 in l|18f) . we have 
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FIG. 4: Dependence of eigenvalues(~ p — 120) on 
anharmonic parameter b. 



non- vanishing matrix elements as 

(4>n, m \pe~ l9 \4>„+i, m+ i) = {<}>n+l,m+l\pe ie \<t>n,m) 

— {n + m + 1) 



\pe 



bn-l,m+l) — {4>n-X,m+l\P& 1 \4>n,m) 



-(n-m-1) 



1/2 



As a result, the matrix elements for (|18f) are 

>) = k(4> nim \pe ze \4> n ^ m >} 



($ + \Hj tI$~ 



= k{— [n±{m- l)]} 1/2 <S„/,„ T i<W,m-i (20) 

ZisJ 



Iff/ T $ / /! 

\*n,ml J-i \ n'm' I 



\pe 



= k{— [n ± (m+ l)]} 1/2 ^,„±i^, m+ i. (21) 

If we assign the quantum numbers j = ±1 to $^ m , Hj-t 
without the anharmonic term connects the states with 
the same quantum number, £ = m — (1/2) j (j = ±1). 
As discussed by Longuet-Higgins|llj. the present matrix 
decomposes into matrices labeled by quantum number 
I, For any given value of £, m can take two values, 
Tji = £—1/2 and £+1/2 corresponding to j = —1 and +1, 
respectively. Thus, the p-th eigenfunction for a given £ 
is expressed as & p ,£- If we consider the trigonal field J^Jl, 
levels with £±3N(N= 1, 2, 3, . . .) are coupled to levels of 
£, as discussed in [Til • The energy matrix is decomposed 
into only three irreducible presentations, A g , B g and E g . 
By the exact diagonalization of each submatrix for the 
Hamiltonian including Hj-t, we can get eigenvalues and 
eigenvectors. The result depends on two parameters, the 
coupling constant k and the strength of the trigonal field 
b. Here, we concentrate on the nearest-neighbor level 
spacing distribution P(s), which plays a prominent role 
in the quantum description of classically chaotic quantum 
systems. The random matrix theory presents a natural 
framework for describing fluctuation properties of spec- 
tra of quantum systems, whose corresponding classical 
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FIG. 5: (a) Histograms of level spacing distribu- 
tion; (b) Integrated level spacing distribution for 
k — 4, b — 5(with trigonal field). Curves of Pois- 
son distribution and Wigner distribution are also 
drawn by open rectangles and crosses, respectively. 



model exhibits chaotic behaviors. In fact the correlations 
in Gaussian ensembles of random matrices are found to 
match very closely the empirical correlations among en- 
ergy levels in classically chaotic systems. If the phase 
space is totally chaotic, the distribution of level spacings 
is Wigner-like(GOE), 



Pw(s) = |sexp (~^s 2 ) • 



By contrast, in classically regular quantum systems the 
levels are independent of each other, and therefore the 
spacings obey Poisson distribution, 

P p (s) = exp(-s). 

Considering these standard criteria, we present the level 
spacing distribution for the Hamiltonian in the fol- 
lowing. 

First, we consider the system without trigonal field(6 = 
0). For the manifold of £ = 1/2 ± 3m(m = 0, 1, 2 . . . , ), 
FigEI a ) shows level spacing distribution characterized by 
the Poisson distribution, which tells that this system is 
regular in the classical dynamics^. Figure 01b) shows 
integrated values of histograms in Figl^a), which well 
fits the function 1 — e~ s . 

Then we reveal the role of trigonal field. Figure0|shows 
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the dependence of eigenvalues on the parameter b, where 
a multitude of avoid crossings can be found. 

Figure Ufa) shows a level spacing distribution in the 
case of b = 5. It should be noted that the distribution 
perfectly agree with the Wigner one with P(s) — * s as 
s — > 0. In Fig|SJb) we show the corresponding integrated 
level spacing distribution, which also well fits the Wigner 
one (GOE). We would like to point out the following fact: 
In the corresponding classical system with k = 4 and 
6 = 5 without use of the adiabatic approximation, the 
almost whole regions in the phase space exhibit stochas- 
tic behavior. However, the system under consideration 
is generic: the phase space of the underlying classical 
dynamics consists of both regular Kolmogorov-Arnold- 
Moser(KAM) tori and chaos. In fact, for smaller values 
of 6(0 < b < 1), we have the Brody distribution interpo- 
lating the Poisson and Wigner limits. 

While most of studies on quantum chaos have been 
limited to the analyses of level-spacing distributions and 
of wavefunctions scars, we shall here embark upon the 
investigation of experimentally accessible new indicators. 
We pay attention to magnetism in the dynamical Jahn- 
Teller system. In particular we discuss the magnetic g- 
factor in the next Sectionl26l. 



given as 



V. ELECTRONIC ORBITAL ANGULAR 
MOMENTUM 



The essential features of the Jahn- Teller coupled sys- 
tems are caused by the interaction between the nuclear 
motions and the electrons in non-Kramers doublets for 
multiplets of d-levels. In particular, the quenching of 
the electronic orbital angular momentum —the so-called 
magnetic g-factor— in those doublets is the fundamen- 
tal subject in magnetism of transition-metal ionic com- 
pounds. As well known, the static Jahn- Teller effect re- 
moves the ground state degeneracy of electronic states. 
This effect leads to the complete disappearance of the 
orbital angular momentum, leaving only spin degree of 
freedom remains. On the other hands the dynamic Jahn- 
Teller systems for the relatively weak coupling have a 
possibility of the nonvanishing orbital angular momen- 
tum in non-Kramers doublets because of the continuous 
distortion of lattices. In fact, in the absence of a trigonal 
field the E g (g> e g system has the freedom of the continu- 
ous distortion of lattices along the continuous minima of 
the adiabatic "Mexican hat" potential. 

Washimiyaj^ payed his attention to the dependence 
of expectation value of this orbital angular momentum 
(L z ) in the excited levels ^ p ,e derived in Section IV. 
He pointed out the oscillatory behavior of (L z ) with in- 
creasing energy levels. The non- vanishing values (L z ) are 



(L z ) p = (*p^=i/ 2 |L»|*p,<=i/a> 



n=l 



(22) 



p=l,2,3,- 



for the vibronic state of I — 1/2. Here, the co- 

efficients of the harmonic oscillator functions. As already 
mentioned, the vibronic wavefunction is given by 

^ P ,i = ai p u-(r)(f> li0 (p,6) + a2 P u + (r)<f)2 < i(p,6) 

+ a 3p u_(r)<j) 3fi (p, 6) + a 4p M+04,i(/5, 0) H (23) 

In (1221 , 5 is the elements in the principal diagonal of the 
matrix. (3 = {u+\L z \u+) = — (u-\L z \u-).) In short, the 
differences of |a ra ,p| 2 between even- and odd- numbers of 
n in the eigenfunction labeled by p play an essential role. 
As p increases, the oscillatory behavior of the angular 
momentum is found. The periods of this oscillation are 
far long in comparison with the variation by odd and even 
numbers for the small coupling k. With increasing k, the 
absolute value of the angular momentum decreases. The 
origin of this oscillation has not been clarified up to now, 
though Washimiya's finding is essential for understanding 
properties of this vibronic systems. 

In what follows, we show a more detailed calculation 
of this oscillatory behavior. Furthermore, the effect of 
the trigonal field is discussed in consideration that this 
field destroys the continuous circular symmetry and that 
the angular momentum £ is not a good quantum num- 
ber. Under that condition, we introduce the density of 
g-factor g{e)de for the electronic orbital angular momen- 
tum in the energy range between e and e + ds as 



g(e)de = ^ 



l t=l/2,n,p 



n=l 



de. 



(24) 



Here, the summation of p is taken over the corresponding 
energy range. It should be noted that levels for £' = 
1/2 ± 3N(N = 1,2,...) are mixed with the levels for 
I = 1/2 in the presence of the trigonal field. 

We show calculated results of g(e) at k = 0.707 in 
FigHj] In FigUJa), the regular oscillatory behavior of 
g(e) for b = as a function of level p, which was re- 
ported by Washimivapt^ three decades ago, is repro- 
duced. In addition to histograms for g(e), we also show 
the envelop-functions constructed by Gaussian coarse- 
graining of each peak. Figures EIb),(c) and (d) show 
the e dependence of g{e) in the presence of trigonal field 
with b = 0.2, 0.3 and 1.41, respectively. Here, we can 
find the suppression of regular oscillation with increasing 
b. 

The emergence of irregular oscillation and the suppres- 
sion of <?-factor with increasing b reflect the underlying 
chaotic behavior in Section III. We find the important 
fact: while the level statistics can reach the Wigner distri- 
bution for sufficiently-large value 6, the suppression of the 
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regular oscillation of g-factor easily occurs for relatively 
small value b. In fact the values of b = 1.41, k = 0.707 
(Q/u) = 1, k x b = 1) in FigEId) are much less than 
the values 6 = 5, k = 4 that guarantee the Wigner distri- 
bution in the case without the adiabatic approximation. 
As a result, the irregular oscillation of g-factor is a pre- 
cursor of quantum chaos, namely, the suppression of its 
regular oscillation occurs even when the classical phase 
space accommodates a partial chaos. 

VI. SUMMARY AND DISCUSSIONS 

We have examined the Jahn- Teller E g tg>e g system from 
a viewpoint of classical chaos and its quantization. Both 
the systems (A) without an trigonal anharmonic term 
and (B) with it are investigated. 

The classical phase space is strictly regular for (A) and 
nonintegrable and chaotic for (B). In general, the system 
(B) is mixed of regular KAM tori and chaos. The rela- 
tive fraction of chaos in the phase space increases as the 
energy or the strength of the anharmonicity is increased. 
In the adiabatic approximation, the full chaos can occur 
easily because of enhancement of non-linear effects due 
to the strong constraint. Without such an approxima- 
tion, we can find the full chaos for relatively large values 



of k and 6, and obtain only the partial chaos, if we adopt 
small values to k and b. 

For the corresponding quantum systems, the level 
spacing distributions are shown to be of Poisson and 
Wigner type for (A) and (B), respectively. We find that 
the Wigner distribution is available as well by adopting 
the relatively large values of k and b. On the other hand, 
the dependence of (/-factors on energy £ is a quite sen- 
sitive indicator of the symptom of chaos in comparison 
with the level spacing distribution. In the system (A) 
g-factor shows regular oscillation with respect to exci- 
tation energy. By contrast, in the system (B) it shows 
a quenched irregular oscillation for the relatively small 
values of k and b. Therefore, we propose the magnetic 
g-factor as a new precursor of quantum chaos, superior 
to the level-spacing distribution. We hope these predic- 
tions will be verified in future such as in the experiment 
of magnetic circular dichromism. 

It will be quite important that the precursor of chaos 
shows up in the observable orbital angular momentum in 
the dynamical Jahn- Teller systems for transition-metal 
ionic-compounds. There are other interesting themes in 
this system, such as an effect of the chaos on a spectral 
feature of phonon side bands, which will also be examined 
in due course. 
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FIG. 6: g-factor g(e)de(e = 0.25) for the electronic orbital angular momentum. £ = 1/2 and k = 0.707. (a),(b),(c) and 
(d) correspond to b — 0,0.2,0.3 and 1.41, respectively. Envelop-functions composed for Gaussian distribution types are also 
depicted. 



